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1-^ I The form factors of weak decays of the B meson to orbitally excited scalar, axial 

Ft ■ vector and tensor light mesons are calculated in the framework of the QCD-motivated 

^ ■ relativistic quark model based on the quasipotential approach. Relativistic effects 

are systematically taken into account. The form factors are expressed trough the 
overlap integrals of the meson wave functions and their dependence on the momen- 
i^ \ tum transfer is selfconsistently determined in the whole kinematical range. On this 

I basis semileptonic and two-body nonleptonic B decay rates to orbitally excited light 



X 



mesons are calculated. Good agreement of the obtained predictions with available 



^ ■ experimental data is found. 

Cs| ■ PACS numbers: 13.20.He, 12.39.Ki 
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INTRODUCTION 



t^^ ■ Recently significant experimental progress has been achieved in studying weak decays of 

O ■ B mesons into hght mesons [L] . Many new decay modes have been measured including non- 

leptonic decays involving excited light mesons |2|-|5|. In Ref. |6| we investigated semileptonic 
B decays to the ground state vr and p mesons in the framework of the relativistic quark model 
based on the quasipotential approach in QCD. The peculiar feature of the heavy-to-light 
weak decays is the very broad kinematical range in which the recoil momentum of the final 
c^ i light meson is significantly larger than the mass of the light meson except the small region 

near the point of zero recoil. Therefore it is very important to take into account all relevant 
relativistic effects and determine the decay form factors without any extrapolations or addi- 
tional parameterizations. Weak decay form factors were selfconsistently calculated |6| in the 
whole accessible kinematical range. It was found that both the behaviour of the form factors 
on momentum transfer and differential and total semileptonic decay rates agree well with 
the rather precise experimental data. This allowed us to determine the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix element Vut- 

In this paper we further extend our analysis for the consideration of the weak semileptonic 
and two-body nonleptonic B decays to the orbitally excited light mesons. We calculate 
the corresponding decay form factors paying special attention for a consistent account of 
all relativistic effects and determination of the form factor dependence on the momentum 
transfer in the whole accessible kinematical range. For calculations we use masses and wave 
functions of orbitally excited light mesons which were previously studied by us in Ref. |7|. 
Table H] quotes a comparison of our predictions for the P-wave light unfiavoured mesons with 



TABLE I: Predicted |3] and measured masses of the P-wave Hght {q = u, d) unflavored mesons (in 

MeV). 
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experimental data [l|. Results for scalar oq, /o, spin triplet i^Pi) axial vector ai, /i, spin 
singlet (^-Pi) axial vector 6i, hi and tensor 02, /2 mesons are presented. Our model predicts 
[TJ the light scalar meson masses heavier than 1 GeV. The scalar mesons below 1 GeV 
are well described as the light tetraquarks, composed from the light scalar diquarks and 
antidiquarks |8[. From Table [I] we see that the calculated mass of the scalar qq state I^Pq 
is significantly lower than the mass of the experimentally observed ao(1450) meson. Our 
model predicts that this state is also the tetraquark composed from the axial vector diquark 
and antidiquark with the mass 1480 MeV [8| . Nevertheless, for the purpose of comparison, 
we here assume that the lightest scalar qq state corresponds to the ao(1450) meson and 
test this assumption by confronting the obtained predictions for the B decays involving this 
meson with available experimental data and results of the different theoretical investigations 
based on quark models, sum rules, light cone sum rules and perturbative QCD. Therefore 
the study of weak B decays involving light scalars is an important problem, since it can help 
to reveal their real nature. 

The calculated weak decay form factors are used for the evaluation of the semileptonic B 
decay branching fractions. Then they are employed for studying the two-body nonleptonic 
B decays using the factorization approach. Such approximation significantly simplifies cal- 
culations, since it expresses the matrix elements of the weak Hamiltonian responsible for the 
nonleptonic decays through the product of the transition matrix elements and meson decay 
constants. Comparison of the obtained results with experimental data, which are mostly 
available for the nonleptonic B decays involving light axial vector mesons |3|, \^ , can help in 
testing this approach and discriminate between different models for form factors. 

II. RELATIVISTIC QUARK MODEL 

In the quasipotential approach a meson is described as a bound quark-antiquark state 
with a wave function satisfying the quasipotential equation of the Schrodinger type 
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where the relativistic reduced mass is 



j7A-sViP,^-,M)^MicD, (1) 
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and El , E2 are the center of mass energies on mass shell given by 

M^ — 1711 + mj ^^ — mf + ml . . 

^' = 2M ' ^' = 2M • (^) 

Here M = Ei + E2 is the meson mass, mi, 2 are the quark masses, and p is their relative 
momentum. In the center of mass system the relative momentum squared on mass shell 

reads 

2 . , [M^ - (mi + m2)^][M^ - (mi - m^)^] 

The kernel ^(p, q; M) in Eq. ([1]) is the quasipotential operator of the quark-antiquark in- 
teraction. It is constructed with the help of the off-mass-shell scattering amplitude, projected 
onto the positive energy states. Constructing the quasipotential of the quark-antiquark in- 
teraction, we have assumed that the effective interaction is the sum of the usual one-gluon 
exchange term with the mixture of long-range vector and scalar linear confining potentials, 
where the vector confining potential contains the Pauli interaction. The quasipotential is 
then defined by [7(] 

V{p, q; M) = Ui{p)u2{-p)V{p, q; M)ui{q)u2{-q), (5) 

with 

V(p, q; M) = ^asD,Ml^l2 + V^U^)TlT2;, + Konf (k), 

where a^ is the QCD coupling constant, D^^, is the gluon propagator in the Coulomb gauge 

Z)°°(k) = -g, D^k) = -^ [s^^ ~ ^y D'^=D^' = 0, (6) 

and k = p — q. Here 7^ and u{p) are the Dirac matrices and spinors 
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where cr and x^ are Pauli matrices and spinors and ej(p) = yp^ + mf . The effective long- 
range vector vertex is given by 

r.(k)=7. + £cr,./^^ (8) 

where k is the Pauli interaction constant characterizing the long-range anomalous chromo- 
magnetic moment of quarks. Vector and scalar confining potentials in the nonrelativistic 
limit reduce to 

Konf(^) = {l-e){Ar + B), 

Kif(r) = e{Ar + B), (9) 

reproducing 

Konf (r) = Konf (0 + Konf (0 = Ar + B , (10) 

where e is the mixing coefficient. 
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FIG. 1: Lowest order vertex function rj( (p, q) contributing to the current matrix element (jlip . 

The expression for the quasipotential of the heavy quarkonia, expanded in f ^/c^ can be 
found in Ref . [9| . The quasipotential for the heavy quark interaction with a light antiquark 
without employing the nonrelativistic (f/c) expansion is given in Refs. [7|, [lO|. All the 
parameters of our model like quark masses, parameters of the linear confining potential 
A and 5, mixing coefficient e and anomalous chromomagnetic quark moment k are fixed 
from the analysis of heavy quarkonium masses and radiative decays [9| . The quark masses 
mb = 4.88 GeV, rric = 1.55 GeV, nig = 0.5 GeV, rriu^d = 0.33 GeV and the parameters of the 
linear potential A = 0.18 GeV^ and B = —0.30 GeV have values inherent for quark models. 
The value of the mixing coefficient of vector and scalar confining potentials e = —1 has been 
determined from the consideration of the heavy quark expansion for the semileptonic B ^ D 
decays [ll| and charmonium radiative decays [9|. Finally, the universal Pauli interaction 
constant k = — 1 has been fixed from the analysis of the fine splitting of heavy quarkonia 
^Pj- states |9| and the heavy quark expansion for semileptonic decays of heavy mesons 
(ll| and baryons [12|. Note that the long-range magnetic contribution to the potential in 
our model is proportional to (1 + k) and thus vanishes for the chosen value of k = — 1 in 
accordance with the fiux tube model. 



III. MATRIX ELEMENTS OF THE ELECTROWEAK CURRENT 

In order to calculate the exclusive semileptonic decay rate of the B meson, it is necessary 
to determine the corresponding matrix element of the weak current between meson states. 
In the quasipotential approach, the matrix element of the weak current J^ = ^7^(1 — 75)^, 
associated with the 6 — )■ g transition, between a B meson with mass Mb and momentum ps 
and a final P-wave light meson F with mass Mp and momentum pp takes the form [l3| 



{FiPF)\J^\BipB)) = /^^^^p,(p)r,(p,q)M/^p,(q), (11) 

where r^(p, q) is the two-particle vertex function and ^E'mpm ^^^ the meson (M = B,F) 
wave functions projected onto the positive energy states of quarks and boosted to the moving 
reference frame with momentum p^^. 

The contributions to r^(p, q) come from Figs. [Hand [H The leading order vertex function 
rL^ (P; q) corresponds to the impulse approximation, while the vertex function T^^\p,q) 
accounts for contributions of the negative-energy states. Note that the form of the relativistic 
corrections emerging from the vertex function F^^^ (p, q) explicitly depends on the Lorentz 
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FIG. 2: Vertex function F^ (p,q) taking the quark interaction into account. Dashed Unes 
correspond to the effective potential V in ([5]). Bold lines denote the negative-energy part of the 
quark propagator. 



structure of the quark-antiquark interaction. In the leading order of the f ^/c^ expansion for 
B and F only r|j^^(p, q) contributes, while r^^^(p,q) contributes at the subleading order. 
The vertex functions look like 



r«(p, q) = u,{p,)^^{l - ^'>)u,{q,){27if5{v, " q?) 



(12) 



and 
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where the superscripts "(1)" and "(2)" correspond to Figs. [T] and [2|, the subscripts q, q, b 
are the quark indices, k = pg — A; k' = q^ + A; A = p^T^ — p^; 
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Here the quark momenta are expressed through relative momenta q and p as follows [13 
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and n*^*^ are three four-vectors given by 



p 



pipj 



"""■<"' = Is- *« + Af(£ + Af) 



E = Jp'^ + M2. 



The wave function of a final P-wave F meson at rest is given by 

*f(p) = HPs+^Pj)ip) = yi^i^Fi^s+^pAp). 



(14) 



where J and M. are the total meson angular momentum and its projection, while S* = 0, 1 
is the total spin. 'il)pf2s+ipj'^{p) is the radial part of the wave function, which has been 
determined by the numerical solution of Eq. ([1]) in [7|, [lO(]. The spin-angular momentum 
part yg-'^ has the following form 

yi^= E(l-^-^i-^2, Sa,+a2\JM){\au ^ ^2|5^i+a2)ri^-'^^-'^^Xi(^i)X2(a2). (15) 

0-10-2 ^ -^ 

Here {jimi, J2^^2\J M) are the Clebsch-Gordan coefficients, Yf^ are spherical harmonics, 
and x(c") (where a = ±1/2) are spin wave functions, 

X(1/2)=Q, x(-l/2)=(°). 

It is important to note that the wave functions entering the weak current matrix element 
flTTj) are not in the rest frame in general. For example, in the B meson rest frame (ps = 0), 
the final meson is moving with the recoil momentum A. The wave function of the moving 
meson \1'fa is connected with the wave function in the rest frame "^fo = ^f by the 
transformation [l3| 

*fa(p) = Dl/\RYjDl/\RYj^Foip), (16) 

where R^ is the Wigner rotation, L^ is the Lorentz boost from the meson rest frame to a 
moving one, and the rotation matrix Dy^{R) in spinor representation is given by 

(J °)<^(i?rj = 5-i(p,)5,(A)5,(p), (17) 

where 
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^ 2mq \ eq{p) 
is the usual Lorentz transformation matrix of the four-spinor 
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IV. FORM FACTORS OF THE SEMILEPTONIC B DECAYS TO THE 
ORBITALLY EXCITED LIGHT MESONS 

The matrix elements of the weak current Jj^ = &7^(1 — 75)^7 for B decays to orbitally 
excited scalar light mesons (S) can be parametrized by two invariant form factors [iJ] 

{S{pF)\qrb\B{pB)) = 0, 

{S{pF)\qr^,b\Bipn)) = Uiq') {p^s+p^f) + /_(?') {p'b - p'f) , (18) 

where four-momentum transfer q = Pb — Pf, Ms is the scalar meson mass. 

The matrix elements of the weak current for B decays to the spin triplet {^Pi) axial 
vector mesons (AV) can be expressed in terms of four invariant form factors [ij] 

{A{pF)\qrb\B{pB)) = {Mb + MA)hvAq'y' + [hvMVB + ^.^Vf]'^ . (19) 

Mb 

{A{pF)Wl,h\B{pB)) = ^^^^'""'"CPBpPFa, (20) 



where Ma and e^ are the mass and polarization vector of the axial vector meson. The 
matrix elements of the weak current for B decays to the spin singlet i^Pi) axial vector 
mesons are obtained from Eqs. (fT9l) by the replacement of the set of form factors hi{q^) by 

The matrix elements of the weak current for B decays to tensor mesons (T) can be 
decomposed in four Lorentz-invariant structures [iJ] 



{np,)\rrmps)) = J^^'^'-'J^Upb.p.. (21) 

{T{pF)\qi>'l!m(PB)) = [Mb + MT)tA,(i-y''^ 

Mb 

+ [tMiqVB + tA,{qVM,^, (22) 

where M^ and e^" are the mass and polarization tensor of the tensor meson. 

The general structure of the current matrix element flTTl) is rather complicated, since 
it is necessary to integrate both with respect to d^p and d^q. We calculate exactly the 
contribution of the leading vertex function F^^^ (p, q) given by Eq. flT2|) to the transition 
matrix element of the weak current (ITT!) using the ^-function. As a result the contribution 
of r^^^ (p, q) to the current matrix element has the usual structure of an overlap integral 
of meson wave functions and can be calculated exactly in the whole kinematical range. 
The calculation of the subleading contribution r^^)(p, q) is significantly more difficult. The 
heavy quark is present in the initial B meson only. Therefore the expansion in its inverse 
powers retains the dependence on the relative momentum in the energy of the final light 
quark. Such dependence does not allow one to perform one of the integrals in the decay 
matrix element fllip using the quasipotential equation. However the final light meson has a 
large recoil momentum (A = p^;' — p^, lAmaxI = {M^ — M|)/(2Mb) ~ 2.5 GeV) almost 
in the whole kinematical range except the small region near g^ = g^^^^ (|A| = 0). This 
also means that the recoil momentum of the final meson is large with respect to the mean 
relative quark momentum |p| in the meson (~ 0.5 GeV). Thus one can neglect |p| com- 
pared to I A| in the final light quark energy eq{p + A) = Jm"^ + (p + A)^, replacing it by 

eq(A) = Jfn^ + A^ in expressions for the T^^\p,q). This replacement removes the relative 
momentum dependence in the energy of the light quark and thus permits to perform one 
of the integrations in the r(f)(p,q) contribution using the quasipotential equation. This 
contribution is relatively small, since it is proportional to the binding energy in the meson. 
To demonstrate this observation, we show in Fig. [3] leading h^ \Q^) and subleading h)^ \q^) 
terms of the form factor hA{q^) as an example. Contributions of such terms to other form 
factors are similar. Therefore application of heavy quark and large recoil energy expansions 
and the extrapolation of the subleading contribution to the small recoil region introduces 
minor errors. Similar calculations for the weak B decays to light ground state mesons were 
made in Ref. [6|]. There it was shown that such extrapolation introduces uncertainties from 
the vicinity of zero recoil q^ = q^^-^ of less than 1%. It is important to emphasize that 
calculating the form factors we consistently take into account all relativistic contributions 
including the boosts of the meson wave functions from the rest reference frame to the moving 
ones, given by Eq. flT6l) . Recently we performed the similar calculation of the weak form 



factors for the transitions of B^ mesons to the orbitally excited mesons [1J|. Since in this 



calculation the spectator charmed quark was treated without the l/rric expansion and the 
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FIG. 3: The form factor h^iq ) of the i? — )• ai weak transition with leading h\ {q ) and subleading 

,(2) 



h^ [q"^) contributions. 



final active quark was considered to be light in the framework of the approach described 
above, we can use for the present calculation the expressions for the decay form factors given 
in Appendix of Ref. IJ] with obvious replacements. In the limits of infinitely heavy quark 
mass and large recoil energy of the final meson, these form factors satisfy all heavy quark 



symmetry relations |l5l . |16 



For numerical evaluations of the form factors we use the quasi pot ential wave functions 
of the B meson and orbitally excited light mesons obtained in [T], [lO|. Our results for the 
masses of these mesons are given in Table [B With the exception of the oq meson they are in 
good agreement with available experimental data [l|. 

In Fig. m we plot form factors of the weak B transitions to the isovector P-wave light 
mesons (scalar ao, axial vector ai and bi, tensor 02). The calculated values of these form 
factors at g^ = 0, several intermediate points and g^ = q'^^.^ = {Mb — MpY are given in 
Table Hn Since we do not distinguish between isovector and isoscalar mesons, the form factors 
for the corresponding weak decays involving isoscalar mesons coincide with the isovector ones 
up to the factor l/v2, which comes from the flavour function of the light neutral meson. 
As it was argued above, the main source of the uncertainties of our form factor calculations 
originates from the subleading terms (IT3l) . We can conservatively estimate the error arising 
from the application of the heavy quark expansion and the extrapolation of the subleading 
contributions to the small recoil region to be less than 1% in the interval g^ = -i- 14 GeV^ 
and to be less than 4% for the interval q^ = 14 GeV^ -^ g^ax- 

These form factors can be approximated with good accuracy by the following expressions: 
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FIG. 4: Form factors of the B decays to the P-wave hght mesons. 



TABLE II: Calculated values of the form factors of the B decays to the P-wave light mesons at 



q^ = 0, several intermediate points and q'^ = 
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= {MB-MFf. 
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0.27 

2.5 0.31 

5 0.35 

7.5 0.39 

10 0.45 

12.5 0.53 

15 0.63 

9max 0-78 



-0.66 
-0.71 
-0.78 
-0.87 
-0.98 
-1.13 
-1.35 
-1.69 
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-1.32 0.04 



0.20 0.17 
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0.11 0.05 



0.74 
0.84 
0.94 
1.06 
1.22 
1.48 



-0.02 2.05 
-0.06 3.48 



-0.04 
-0.04 
-0.04 
-0.04 
-0.03 
-0.02 
0.01 
0.08 



-0.04 
-0.04 
-0.04 
-0.04 
-0.03 
-0.03 
-0.02 
-0.02 



-0.05 
-0.08 
-0.09 
-0.10 
-0.11 
-0.13 
-0.16 
-0.27 



-0.29 
-0.31 
-0.32 
-0.35 
-0.38 
-0.44 
-0.60 
-1.02 



-0.29 
-0.30 
-0.33 
-0.38 
-0.45 
-0.56 
-0.82 
-1.13 



-0.16 
-0.18 
-0.20 
-0.22 
-0.24 
-0.26 
-0.29 
-0.35 



-0.01 
-0.02 
-0.02 
-0.03 
-0.04 
-0.06 
-0.18 
-0.43 



-0.01 
-0.01 
-0.01 
-0.01 
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where the values -F(O) are given in Table HTl and the values cxj {i = 1,2,3,4) are given in 
Table lllli ^ In the next sections we apply the obtained form factors for the calculation of 
semileptonic and nonleptonic B decays involving orbitally excited light mesons. 



^ Note that we do not use these parameterizations for further calculations. 
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TABLE III: Fitted parameters of form factor parameterizations (j23p and (|24p . 
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0-1 -1.20 -0.63 -0.91 0.04 -0.15 -0.53 1.20 -0.95 -3.04 0.72 0.98 -0.89 -4.25 -1.33 

0-2 -0.27 -2.29 -2.40 -1.59 -1.85 -7.10 -12.5 9.71 2.46 -14.3 -14.7 -4.04 7.39 16.9 

(73 2.35 6.03 9.36 0.36 -0.82 23.9 38.1 -30.9 13.5 45.5 36.8 17.0 -3.10 -63.9 

(74 -2.53 -5.34 -9.90 -0.93 2.52 -24.2 -41.7 37.1 -22.2 -46.4 -31.8 -17.4 -3.69 77.1 



V. SEMILEPTONIC B DECAYS TO ORBITALLY EXCITED LIGHT MESONS 



The differential decay rate for the B meson decay to P-wave hght mesons reads 17 
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where Gp is the Fermi constant, Vub is the CKM matrix element, A = X{M'^,Mp,q'^) = 
M^ + M^ + q^ - 2(M|M| + M^q^ + M^q^), mi is the lepton mass and 

HH^ = H+hI + H_h1_ + HoHI (26) 

The helicity components if-t, Hq and Ht of the hadronic tensor are expressed through the 
invariant form factors. 
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FIG. 5: Predictions for the differential decay rates (in |T/ubplO ^^ GeV ^) of the B semileptonic 
decays to the P-wave hght mesons. 
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Here the subscripts ±,0,t denote transverse, longitudinal and time helicity components, 
respectively. 

Now we substitute the weak decay form factors calculated in the previous section in the 
above expressions for decay rates. The resulting differential distributions for the B decays 
to the P-wave light mesons are plotted in Fig. [5l The corresponding total decay rates are 
obtained by integrating the differential decay rates over q^. For calculations we use the 
following value of the CKM matrix element \Vub\ = 0.0038 ± 0.00044 (l|. It is necessary 
to point out that the kinematical range accessible in these semileptonic decays is rather 
broad. Therefore the knowledge of the g^ dependence of the form factors is very important 
for reducing theoretical uncertainties of the decay rates. Our results for the semileptonic 
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TABLE IV: Predictions for the decay rates of the semileptonic B decays to the P-wave hght mesons 

(in iKbpps-^). 

Decay T{l = e,fi) T{1 = t) 
B -^ aolu 3.18 1.00 

B -^ ailv 5.90 1.95 

B -^ hilv 4.04 0.91 

B -^ a2li' 5.77 1.27 



TABLE V: Comparison of theoretical predictions for the branching ratios of the semileptonic B 
decays to the P-wave light mesons (in 10~^). EFG: Ebert, Faustov, Galkin, this paper. 

Decay EFG ISGW2 [j^ pQCD [19] LCSR [20] SR ^ 

7^3 

19.2 29.6+l^-S 30.2+|!^1 16 

24.1 

1.1 



B decay rates to the P-wave light mesons are given in Table IIVI ^ The errors of the decay 
rates presented in this table originate from the form factor uncertainties, discussed in the 
previous section, and are about 7%. 

The predictions for the branching ratios of the semileptonic B decays to the P-wave 



po -^ a^ev 


7.0 ±2.8 


P° -^ a^TV 


2.2 ±0.9 


P° -^ afeu 


13.0 ±5.2 


po -^ a+Tu 


4.3 ±1.7 


po -^ btev 


8.9 ±3.6 


po ^ btrv 


2.0 ±0.8 


po -^ ateu 


12.7 ±5.1 


P° -^ at TV 


2.8 ±1.1 



09 c;+23.6 
•3 '^•'-'-13.6 


istl 


13.2t|? 


6.3l|t 


29.6+13'i 


3o.2tl°:i 


l^Atll 




no Q+15.1 
ZO.(5_]^2.2 


19.311^ 


12.6+f-^ 




ii.6t«:^ 


16 


^■Itl-l 


6 



light meson states ^ are compared with the previous calculations 18l-l2ll] in Table |Vl The 
authors of Refs. [l8| use the constituent quark (ISGW2) model. Perturbative QCD (pQCD) 
approach is adopted in Ref. (l9|]. Calculations in Ref. [20| are based on the light-cone 
QCD sum rules (LCSR), while Ref. [21!] employs the QCD sum rules. Our predictions for 
B — >■ ao,i(6i)/z/ decays are somewhat lower than the central values of pQCD and LCSR, but 
they are consistent within rather large errors of the considered approaches. The predicted 
central values of the branching ratio for the B — )■ a2li' decay are close in all calculations 
except the ISGW2 model, which predicts an order of magnitude lower value. Therefore we 
find that essentially different theoretical approaches give the values for the B — )■ aj{hi)li> 
decay branching ratios of order of 10"'^ which is the same as for the decays to the ground 
state TT and p mesons. It is important to verify these predictions experimentally. 



^ We used the calculated value of the ao mass from Table HI 

■^ The presented errors of our calculations arise from the theoretical uncertainties in form factor calculations 

and experimental uncertainties which mainly originate from the rather poor knowledge of the CKM matrix 

element Vub- The latter uncertainty is dominant. 
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VI. NONLEPTONIC DECAYS 

In the standard model nonleptonic B decays are described by the effective Hamiltonian, 
obtained by integrating out the heavy H^-boson and top quark. For AB = 1 transitions 
{q = d,s) Jin 

(30) 
The Wilson coefficients Ci^jj) are evaluated perturbatively at the W scale and then are 
evolved down to the renormalization scale /i ~ m^ by the renormalization-group equations. 
The expressions Oi are local four-quark operators which are given by 

Oi = {q'b)v-A{qq')v-A , 

O2 = {q'ibj)v-Aiqjq'i)v-A , 

03(5) = {qb)v-AY.(^'l')v-A,{V+A) , 

<?' 

^4(6) = {qibj)v-A'Yiqjq'i)v-A,iV+A) , 

g' 
3 

^7(9) = -iqb)v-Aj2'^g'(^''i')v+A,{v-A) , 

<?' 
3 

08(10) = -{qibj)v-AY.^g'(^j^'i^V+A,{V-A) , (31) 

g' 

where e^ denotes the quark electric charge and the following notations are used 

{qq')vTA = qit^{'^Ti5)q'- 

The nonleptonic two-body decay amplitude of a i? meson into light mesons can be ex- 
pressed through the matrix element of the effective weak Hamiltonian i^es ^ 

M{B ^ M1M2) = {M^M2\H,s\B) = ^Jk^K; [ci{MiM2\0^\B) + C2{MiM2\0^\B)] 

10 ^ 

-VtbV,lY.('^(M,M2mfi)\B) . (32) 

i=3 J 

The factorization approach, which is extensively used for the calculation of two-body non- 
leptonic decays assumes that the nonleptonic decay amphtude reduces to the product of a 



meson transition matrix element and a decay constant [2J]. This assumption in general can- 
not be exact. However, it is expected that factorization can hold for energetic decays, where 
both final mesons are light and therefore possess large recoil momenta^ 25|. A justification of 



this assumption is usually based on the issue of color transparency 



. In these decays the 



^ Since we make all further calculation adopting the naive factorization assumption we neglect contributions 
of charming penguins [23j which are absent in this approximation. 
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final hadrons, which have a large relative momentum, are produced in the form of almost 
point-like color-singlet objects that do not couple to soft gluons at leading order. Therefore 
they undergo only hard interactions with the B meson remnants before they hadronize. A 



more general treatment of factorization is given in Refs. [27[ where it is shown that its naive 
form follows in the heavy quark limit at zeroth order in as and Aqcd/^^b- 

Then the decay amplitude can be approximated by the product of one-particle matrix 
elements, e.g., the tree part of the matrix element {q = d, s) is given by 

{M^M^\c,0, + C202\B-) ^ a^{M^,\{ub)v-A\B-){M^\iqu)v-A\0) 

+a2{M^\{qb)v-A\B-){M',\{uu)v^A\0), (33) 

in which the Wilson coefficients appear in the following linear combinations 

tti = Ci + -Tj-Ci+i {i = odd) 

at = Ci + iTj-Ci-i {i = even) (34) 

and Nc is the number of colors. For calculations we use the values of the next-to- leading order 



Wilson coefficients obtained in Ref. [22[ for A^^ = 225 MeV in the HV ('t HooftVeltman) 



scheme: Ci = 1.105, Cs = -0.228, C3 = 0.013, C4 = -0.029, Cg = 0.009, Cg = -0.033, 
cj/a = 0.005, Cg/a = 0.060, Cg/a = —1.283, Cio/a = 0.266, where a is the fine structure 
constant. 

The matrix elements of the weak current {M\{qb)v-A\B) between B and light meson 
states are expressed through the decay form factors (see, e.g., flT8|) - fl22|) ). The matrix element 
(M|(gig2)y-yi|0) between vacuum and a final pseudoscalar (P), vector {V), scalar (S) and 
axial vector {AV) meson is parametrized by the decay constants fpy,s,Av 

(Plqil^l^MO) = ^fpp'p, 
{V\q,^,q2\0) = e,Mvfv, 
{S\qirq2\0) = fsp'p, 

(v4\/|gi7^75g2|0) = e^MAvfAv, (35) 

while the corresponding matrix element for the tensor (T) meson vanishes since 

(T|gi7^g2|0)oce^,p'' = 0, (36) 

due to the properties of the polarization tensor e^^. The pseudoscalar fp and vector fv 
decay constants were calculated within our model in Ref. ^]. It was shown that the 
complete account of relativistic effects is necessary to get agreement with experiment for 
decay constants especially of light mesons. The scalar decay constant is proportional to the 
difference of the light quark masses and thus exactly vanishes for the neutral scalar mesons 
and also vanishes for the charged ones if isospin symmetry is assumed, f^± ~ /^o = //^ = 0. 
The axial vector bi (^-Pi) meson cannot be produced from the vacuum by the axial vector 
current due to the G-parity conservation and, therefore, fb^ = 0. We use the following 
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values of the decay constants: U = 0.131 GeV, fp = 0.220 GeV, fx = 0.160 GeV, fx* = 
0.220 GeV, /^ = 0.240 GeV and fa, = ff, = 0.238 GeV. The relevant CKM matrix elements 
Q are |Kd| = 0.975, |Ks| = 0.225, iKfel = |0.0019 - i0.0033| = 0.0038, \Vtb\ = 0.999, 
\Vtd\ = 10.0080 -i 0.00321 = 0.0086, |V"i,| = 0.0403. 

The matrix elements of the weak current between the B meson and the final light meson 
entering the factorized nonleptonic decay amplitude (l33l) are parametrized by the set of 
decay form factors. Using the form factors obtained in Sec. [IV] and our previous results for 
B decays to the ground-state vr and p mesons [6| , we get predictions for the branching ratios 
of the nonleptonic B decay to orbitally excited light mesons and display them in Tables 
IIXI in comparison with other calculations 29l-l40p and available experimental data [l- 



We can roughly estimate the error of our calculations within adopted naive factorization 
approach to be about 40%. It originates both from theoretical uncertainties in the form 
factor and effective Wilson coefficient calculations and experimental uncertainties in the 
values of CKM matrix elements (which are dominant), decay constants and meson masses. 
The difference between charge combinations of the same final mesons in Tables IVItilXI 
is sometimes enormous. It originates mainly from the different set of diagrams for the 
nonleptonic decay process involving charged and neutral light mesons. For example, decays 
B^ — i- aQ2K^ and B^ -^ aj] g-ft'" are tree dominated, while decays B^ — t- Og 2-?^° and 
B^ — 7- aQ2K^ are penguin dominated. In the case of the decays involving the axial vector 
Ci meson the situation is more complicated since its decay constant is not equal to zero 
and thus additional diagrams, where ai is produced by the weak current from the vacuum, 
contribute to the nonleptonic decay amplitudes. It is also necessary to take into account 
additional factors arising form the composition of neutral light unflavored mesons [e.g.. 



a° = 1/V2{dd -uu),i = 0, 1, 2]. 



In Table I VI I we present predictions for the branching ratios of the two-body nonleptonic 
B decays involving the P-wave light scalar qq meson oq. We compare our results with 



predictions of LCSR [29[, QCD factorization with the form factors evaluated in the light- 
front quark mod el [30l| . naive factorization with form factors obtained using QCD sum rules 
3l| . and pQCD [32] approach. We see that LCSR 29| and pQCD |32[ give the branching 



fractions which are almost an order of magnitude larger than our predictions. This is the 
consequence of significantly larger values of the form factors /+"°(0) = 0.46 in LCSR [29| 
and /f'"'(0) = 0.86 in pQCD Q approaches compared to our result /+''°fO) = 0.27. On 
the other hand, our predictions are consistent with the ones of [30] and |3l| which both 
use /+"°(0) = 0.26. In Ref. [30| it was argued that the two-body nonleptonic B decay 
rates involving light scalars are significantly different in two- and four-quark pictures of 
these mesons. Therefore experimental measurement of these nonleptonic decay rates can 
discriminate between different models for form factors and help to clarify the nature of 
light scalars. At present only experimental upper limits [2|] are available for two B^ decay 
modes to the scalar ao(1450) and charged pion or kaon, but unfortunately they involve the 
unmeasured branching ratio i3(ao(1450) — )■ rjir). 

In Tables IVIII and IVIIII we compare theoretical predictions for the two-body nonleptonic 
B decays involving the axial vector ai and 61 light mesons with available experimental data. 
Naive factorization hypothesis and decay form factors calculated within the ISGW2 model 



are used in Ref. [33|. In Ref. [34| these decays are investigated in the framework of QCD 



factorization with the light cone distribution amplitudes evaluated using QCD sum rules. 



Only central values are given for all theoretical predictions which have rather large error bars. 
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TABLE VI: The branching ratios of the two-body nonleptonic B decays involving the scalar I^Pq 
light mesons (in 10~^ ). 

Decay EFG [29] [30^ ^1^ [2| Exp. [1, 2j 

20 3.1 8 < 2.3/B{ao -^ r]7r) 



B^ -^ a^TT- 


3.6 


B^ ^ agTrO 


0.03 


B- -^ agvr^ 


2.0 


B- -^ a^TT^ 


0.07 


B^ -^ a^K- 


0.29 


^0 ^ aC^o 


0.02 


B -^ a^K^ 


0.04 


B- -^ agK- 


0.15 


B^ -^ a+p- 


10.3 


^0 ^ ay 


0.05 


B- ^ ay 


5.5 


B- -^ a^p^ 


0.10 


B° -^ a-^K*- 


2.0 


B^ -^ a[]K*o 


0.35 


B- -^ agK*- 


1.1 


B- -^ a^K*^ 


0.8 


B° -^ a^a{ 


11.4 


B^ ^ ay 


0.05 


B- -^ a^al 


0.11 


B- -^ ay 


6.1 



38 



3.1 


8 


0.7 




2.5 


4 


1.1 


0.01 


0.3 


1 


0.1 




0.1 




0.2 


0.5 


13.3 




3.2 




25.4 




4.5 




5.3 




2.7 




2.6 




7.8 





< 3.l/B{ao -^ rjTr) 



28 
14 
7.0 
30 



The authors of Ref. [35| employ naive factorization and additional input of a limited number 
of experimental data ^, while pQCD and soft collinear effective theory (SCET), with form 
factors being fitted parameters, are applied in Ref. |36j . The two-body nonleptonic B decays 
involving axial vector light mesons are the best studied experimentally among the decays to 
excited light mesons. Values or upper limits are available for almost a half of the decays given 
in Tables I VI It IVIIII Notwithstanding rather large experimental error bars the existing data 
can already discriminate between various theoretical approaches. As it is seen from these 
tables our results and predictions of Refs. [3^, J35|, |36[ (pQCD and SCET) are consistent 
with each other (taking into account rather large error bars) for the most decay branching 
ratios and agree with the available experimental data. While the results of Ref. |33| are in 
most cases significantly different and seem to be ruled out by experiment. 

In Table HXl our predictions for the branching ratios of the two-body nonleptonic B decays 
involving tensor 02 and /2 light mesons are confronted with other theoretical predictions and 



experimental data. References [37[ and |38[ employ generalized factorization complemented 
by form factors calculated in the nonrelativistic ISGW model and covariant light-front ap- 
proach, respectively. The authors of Ref. (39| apply the QCD factorization, while naive 



^ Two possible values of the i4ri(f270) and ii'i(1400) mixing angle are considered in Ref. [35|. In Tables IVTll 
and IVIIll we present results only for its preferred value 58°. 
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TABLE VII: The branching ratios of the two-body nonleptonic B decays involving the axial vector 
I'^Pi light mesons (in 10~^ ). 



Decay 


EFG 


[33] 


[34] 


[35] 


pQCD [36] 


[36] 


Exp. [1, 3] 


5° -^ afTT^ 


15.7 


74.3 


9.1 


11.8 


12.7 


10.7 


13.0 ±4.3 


B° -^ af 7r+ 


21.1 


36.7 


23.4 


12.3 


15.7 


17.0 


24.2 ±5.8 


B^ ^ a±7r± 


36.8 


111.0 


32.5 


24.1 


28.3 


27.7 


33 ±5 


B° ^ a?7r0 


0.34 


0.27 


0.9 


1.7 


0.12 


5.5 


< 1100 


B- -^ a'j'vr" 


11.3 


43.2 


7.6 


8.8 


6.7 


17.2 


20.4 ±4.7 ±3.4 


B- -^ aj^vr" 


13.7 


13.6 


14.4 


10.6 


8.1 


19.0 


13.2 ±2.7 ±2.1 


B° -^ afK- 


13.2 


72.2 


18.3 


41 


20.6 


15.8 


16.3 ±2.9 ±2.3 


^0 ^ a?KO 


6.0 


42.3 


6.9 


25 


8.0 


6.3 




5- ^ af^o 


19.8 


84.1 


21.6 


52 


25.5 


15.5 


33.2 ±5.0 ±4.4 


B- -^ a\K~ 


14.3 


43.4 


13.9 


28 


15.4 


10.5 




B° -^ a^p- 


20.7 


4.3 


23.9 








<61 


B^ -^ aiK*- 


3.9 


0.92 


10.6 










B- -^ a^K*^ 


0.66 


0.51 


11.2 








-, q+l. 1 + 1.1 


PO ^ a?0 


0.001 


0.0005 


0.01 










B° -^ a^a{ 


46.1 


6.4 


37.4 








47.3 ±10.5 ±6.3 


B^ ^ a^,a^^ 


0.81 


0.1 


0.5 










B~ — ;> afa? 


31.5 


3.6 


22.4 








< 13000 


^° ^ a?/i 


0.85 


0.02 


0.1 










B- ^ a^/i 


17.9 


3.7 


12.4 










B' ^ /ivrO 


0.56 


0.47 


0.26 










S- ^ /ivr- 


11.6 


34.1 


5.2 










5° ^ /li^o 


2.9 


34.7 


14.6 










5- ^ /li^- 


4.9 


31.1 


14.8 








<2.0 



factorization and the improved ISGW2 model are used in Ref. [40[. From this table we 
see that theoretical predictions strongly depend on the adopted approach and the model 
for form factors. Experimental data are available only for a few considered decay modes 
and represent mostly upper limits. The measurements [5J] were recently carried out for two 
charged B decays involving the tensor /2 meson and the charged pion and kaon as well as for 
one neutral B decay to /2 and K^. Our model prediction for B^ — t- f2K~ is in agreement 
with experiment, while the ones for B^ — > /2-ft'° and for B^ -^ /27r~ are lower and larger 
than experimental values, respectively. However experimental errors are still large in order 
to make definite conclusions. 



VII. CONCLUSIONS 

Weak form factors of the B mesons decays to the first orbital excitations of light mesons 
were calculated in the framework of the QCD-motivated relativistic quark model. The form 
factor dependence on the momentum transfer was selfconsistently determined in the whole 
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TABLE VIII: The branching ratios of the two-body nonleptonic B decays involving the axial vector 
l^Pi light mesons (in 10^^ ). All experimental values include the unmeasured branching ratios 



Decay 


EFG 


[33] 


m 


[35] 


pQCD [36] 


[36] 


Exp. [1, 4] 


B° -^ b+TT- 


17.7 


36.2 


11.2 


0.7 


18.7 


7.7 




B^ -^ b{TT+ 








0.3 


«0 


1.4 


0.6 




B^ ^ bfn^ 


17.7 


36.2 


11.5 


0.7 


20.2 


8.3 


10.9 ±1.2 ±0.9 


B^ -^ 5?7r° 


0.18 


0.15 


1.1 


0.01 


1.5 


1.8 


0.4 ±0.8 ±0.2 


B- -^ b^TT- 


9.5 


18.6 


9.6 


0.7 


5.1 


5.0 


6.7 ±1.7 ±1.0 


B -^ b{7r° 


0.62 


0.29 


0.4 


0.5 


1.0 


2.0 


1.8 ±0.9 ±0.2 


B° -^ b^K- 


11.6 


35.7 


12.1 


2.0 


42.9 


8.5 


7.4 ±1.0 ±1.0 


B^ ^ b^^K^ 


4.4 


19.3 


7.3 


4.0 


23.3 


4.0 


5.1 ±1.8 ±0.5 


B- -^ b{K° 


8.3 


41.5 


14.0 


3.0 


55.0 


8.6 


9.6 ±1.7 ±0.9 


B- -^ blK- 


8.2 


18.1 


6.2 


0.7 


24.9 


4.6 


9.1±1.7±1.0 


B' ^ btp- 


22.7 


1.6 


32.1 










50 ^ 50pO 


0.20 


0.002 


3.2 








<3.4 


B- ^ 6rp° 


0.43 


0.0005 


0.9 








<5.2 


B- -^ b\p- 


11.2 


0.86 


29.1 








<3.3 


B° -^ bfK*- 


4.1 


0.32 


12.5 










BO ^ 50^*0 


0.36 


0.15 


6.4 








<8.0 


B- -^ blK*- 


2.2 


0.12 


12.8 








<6.7 


B- -^ b^K*° 


0.7 


0.18 


7.0 








<5.9 


50 ^ 5?(/> 


0.001 


0.0002 


0.01 










50 -^ b+b^ 








1.0 










^0 ^ b\b\ 








3.2 










B- -^ b{b1 








1.4 










B^ ^ /ii^o 


0.22 


0.16 


0.16 










B^ — ^ hiTT^ 


9.6 


18.6 


1.8 










B^ ^ hiK^ 


4.3 


19.0 


10.9 










B- -^ hiK- 


11.1 


19.0 


11.3 











accessible kineinatical range without applying any additional paranieterizations and extrap- 
olations. All relativistic contributions, including contributions of the intermediate negative- 
energy states and transformations of the wave functions to the moving reference frame were 
consistently taken into account. This significantly reduces theoretical uncertainties of the 
obtained form factors. 

On this basis the branching ratios of the B semileptonic decays to orbitally excited light 
mesons were calculated. Our predictions were compared with other theoretical calculations 
based on the ISGW2 quark model [l8|, perturbative QCD [l9J, light cone sum rules |20| and 



QCD sum rules [21|. It is important to point out that in most of the previous approaches 
the weak form factors were calculated in some particular kinematical point or limited kine- 
matical range and than were extrapolated to the whole accessible kinematical range, which 
is rather broad for such decays. Thus the ISGW2 quark model allows the calculation of 
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TABLE IX: The branching ratios of the two-body nonleptonic B decays involving the tensor 1^P2 
Hght mesons (in 10~^ ). 



Decay 


EFG 


m 


[38] 


[39] 


[40] 


Exp. [1, 5] 


B^ -^ a^TT- 


9.8 


4.9 


8.19 


5.2 


13.0 


<300 


B^ ^ aOvrO 


0.009 


0.0003 


0.007 


0.24 


0.18 




B- -^ a^vr" 


5.2 


2.6 


4.38 


3.0 


6.7 




B- -^ a^7r° 


0.19 


0.001 


0.015 


0.24 


0.38 




5° -^ a^K- 


1.6 


0.58 


0.73 


9.7 


0.95 




^0 ^ a^K" 


0.02 


0.005 


0.014 


4.2 






B- -^ a^K^ 


0.05 


0.011 


0.015 


8.4 






B -^ a^K- 


1.0 


0.31 


0.39 


4.9 


0.51 


<45 


B^ -^ a^p- 


27.1 


14.7 


36.2 


11.3 


36.2 




B^ ^ ay 


0.23 


0.003 


0.03 


0.39 


0.5 




B- -^ a^p- 


14.6 


7.3 


19.3 


8.4 


19.4 




B- -^ a^p^ 


0.51 


0.007 


0.071 


0.82 


1.1 


<720 


B° -^ a^K*- 


5.0 


3.5 


7.25 


6.1 


1.9 




BO ^ aO;^*o 


0.9 


2.1 


4.0 


3.4 






B- -^ a^K*^ 


1.8 


4.5 


8.6 


6.1 






B- -^ 4k*- 


2.7 


1.9 


2.8 


2.9 


1.0 




B^ -^ a^a{ 


59.0 




42.5 








B^ ^ ay^ 


0.32 




0.04 








B —?- 020.1 


0.69 




0.085 








B —?- O2O1 


31.6 




22.7 








B^ ^ /2^° 


0.09 


0.0003 




0.15 


0.19 




B ^ f2TT- 


5.1 


2.8 




2.7 


7.1 


1.57 ± 0.42 ±0.16l|j:f^ 


50 ^ /2i^0 


0.11 


0.005 




3.4 




2.7tJ|±0.9 


B- -^ f2K- 


1.0 


0.34 




3.8 


0.54 


1.33 ± 0.30 ±0.1ll|jJ^ 



the form factors at g^ = q"^^^ and then appUes the Gaussian parameterization for them, 
while light cone sum rules determine form factors in the range near g^ = and, therefore, 
require extrapolation. It was found that all these essentially different approaches predict 
that semileptonic decays to orbitally excited light mesons have branching ratios of order 
10~^ which is the same as for the decays to the ground state vr and p mesons. 

The obtained form factors were used for the evaluation of the branching ratios of the 
two-body nonleptonic decays of B mesons involving orbitally excited light mesons. The 
factorization approach was employed to express the decay matrix elements through the 
products of the weak form factors and decay constants. Decays involving scalar ao(1450),^ 
axial vector ai(1260), /i(1285), 6i(1235), /ii(1170) or tensor a2(1320), /2(1270) and light vr, 
p, K and K* mesons were considered. Obtained predictions were compared with previous 
calculations based on naive and generalized factorization with form factors obtained in dif- 



^ Under the assumption that it is the I^Pq qq state. 
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ferent models, QCD factorization, light cone sum rules and perturbative QCD. It was found 
that the results significantly depend on the adopted approach for the calculation of the decay 
matrix elements and form factors. Our predictions agree well with the experimental data, 
which are mostly available for the decays involving axial vector light mesons, while some of 
the previous calculations significantly deviate from experimental values. Future more pre- 
cise and comprehensive data, especially on the semileptonic decays, can help to discriminate 
between various theoretical approaches and form factor models. 



Acknowledgment s 

The authors are grateful to M. Miiller-Preussker for support and to V. Matveev, V. Savrin 
and M. Wagner for discussions. Two of us (R.N.F. and V.O.G.) acknowledge the support 
by the Deutsche Forschungsgemeinschaft under contract Eb 139/6-1. 



[1] K. Nakamura [Particle Data Group], J. Phys. G 37, 075021 (2010). 
[2] B. Aubert et al. [BABAR Collaboration], Phys. Rev. D 75, 111102 (2007). 
[3] B. Aubert et al. [BABAR Collaboration], Phys. Rev. Lett. 98, 181803 (2007); Phys. Rev. Lett. 
99, 261801 (2007); Phys. Rev. Lett. 100, 051803 (2008); Phys. Rev. D 80, 092007 (2009); 
P. del Amo Sanchez et al. [The BABAR Collaboration], Phys. Rev. D 82, 091101 (2010). 
[4] B. Aubert et al. [The BABAR Collaboration], Phys. Rev. Lett. 99, 241803 (2007); Phys. Rev. 

D 78, 011104 (2008); Phys. Rev. D 80, 051101 (2009). 
[5] A. Garmash et al. [Belle Collaboration], Phys. Rev. Lett. 96, 251803 (2006); B. Aubert et al. 
[BABAR Collaboration], Phys. Rev. D 78, 012004 (2008); Phys. Rev. D 79, 072006 (2009); 
Phys. Rev. D 80, 112001 (2009). 
[6] D. Ebert, R. N. Faustov and V. O. Calkin, Phys. Rev. D 75, 074008 (2007). 
[7] D. Ebert, R. N. Faustov and V. O. Calkin, Phys. Rev. D 79, 114029 (2009). 
[8] D. Ebert, R. N. Faustov and V. O. Calkin, Eur. Phys. J. C 60, 273 (2009). 
[9] D. Ebert, R. N. Faustov and V. O. Calkin, Phys. Rev. D 67, 014027 (2003). 
[10] D. Ebert, V. O. Calkin and R. N. Faustov, Phys. Rev. D 57, 5663 (1998) [Erratum-ibid. D 59, 

019902 (1999)]; D. Ebert, R. N. Faustov and V. O. Calkin, Eur. Phys. J. C 66, 197 (2010). 
[11] R. N. Faustov and V. O. Calkin, Z. Phys. C 66, 119 (1995). 
[12] D. Ebert, R. N. Faustov and V. O. Calkin, Phys. Rev. D 73, 094002 (2006). 
[13] R. N. Faustov, Ann. Phys. 78, 176 (1973); Nuovo Cimento A 69, 37 (1970). 
[14] D. Ebert, R. N. Faustov and V. O. Calkin, Phys. Rev. D 82, 034019 (2010). 
[15] J. Charles, A. Le Yaouanc, L. Oliver, O. Pene and J. C. Raynal, Phys. Rev. D 60, 014001 

(1999). 
[16] D. Ebert, R. N. Faustov and V. O. Calkin, Phys. Rev. D 64, 094022 (2001). 
[17] M. A. Ivanov, J. C. Korner and P. SantorelH, Phys. Rev. D 71, 094006 (2005) [Erratum-ibid. 

D 75, 019901 (2007)]. 
[18] D. Scora and N. Isgur, Phys. Rev. D 52, 2783 (1995). 

[19] R. H. Li, C. D. Lu, W. Wang and X. X. Wang, Phys. Rev. D 79, 014013 (2009); R. H. Li, 
C. D. Lu and W. Wang, Phys. Rev. D 79, 034014 (2009); W. Wang, Phys. Rev. D 83, 014008 
(2011). 



21 



[20] Y. M. Wang, M. J. Aslam and C. D. Lu, Phys. Rev. D 78, 014006 (2008); K. C. Yang, Phys. 



[21 
[22 
[23 
[24 
[25 
[26 
[27; 

[28 
[29 
[30 

[31 

[32 
[33 

[34; 

[35; 

[36 

[37; 

[38 
[39; 

[4o; 



Rev. D 78, 034018 (2008); Z. G. Wang, arXiv:1011.3200 [hep-ph]. 

T. M. Aliev and M. Savci, Phys. Lett. B 456, 256 (1999). 

G. Buchalla, A. J. Buras and M. E. Lautenbacher, Rev. Mod. Phys. 68, 1125 (1996). 

M. Ciuchini, E. Franco, G. Martinelh and L. Silvestrini, Nucl. Phys. B 501, 271 (1997). 

M. Bauer, B. Stech, and M. Wirbel, Z. Phys. C 34, 103 (1987). 

M. J. Dugan and B. Grinstein, Phys. Lett. B 255, 583 (1991). 

J. D. Bjorken, Nucl. Phys. B (Proc. Suppl.) 11, 325 (1989). 

M. Beneke, G. Buchalla, M. Neubert and C. T. Sachrajda, Phys. Rev. Lett. 83, 1914 (1999); 

Nucl. Phys. B 591, 313 (2000). 

D. Ebert, R. N. Faustov and V. O. Galkin, Phys. Lett. B 635, 93 (2006). 

V. Chernyak, Phys. Lett. B 509, 273 (2001). 

H. Y. Cheng, C. K. Chua and K. C. Yang, Phys. Rev. D 73, 014017 (2006); Phys. Rev. D 77, 

014034 (2008). 

D. Delepine, J. L. Lucio M., J. A. Mendoza S. and C. A. Ramirez, Phys. Rev. D 78, 114016 

(2008). 

Z. Q. Zhang, Phys. Rev. D 83, 054001 (2011). 

G. Calderon, J. H. Munoz and C. E. Vera, Phys. Rev. D 76, 094019 (2007). 

H. Y. Cheng and K. C. Yang, Phys. Rev. D 76, 114020 (2007); H. Y. Cheng and K. C. Yang, 

Phys. Rev. D 78, 094001 (2008) [Erratum-ibid. D 79, 039903 (2009)]. 

V. Laporta, G. NardulH and T. N. Pham, Phys. Rev. D 74, 054035 (2006) [Erratum-ibid. D 

76, 079903 (2007)] 

W. Wang, R. H. Li and C. D. Lu, Phys. Rev. D 78, 074009 (2008). 

C. S. Kim, J. P. Lee and S. Oh, Eur. Phys. J. C 22, 683 (2002); Eur. Phys. J. C 22, 695 

(2002); Phys. Rev. D 67, 014002 (2003). 

J. H. Munoz and N. Quintero, J. Phys. G 36, 095004 (2009); J. Phys. G 36, 125002 (2009). 

H. Y. Cheng and K. C. Yang, Phys. Rev. D 83, 034001 (2011). 

N. Sharma, R. Dhir and R. C. Verma, Phys. Rev. D 83, 014007 (2011). 



